Quantum Simulation of the Ultrastrong Coupling Dynamics in Circuit QED 
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We propose a method to get experimental access to the physics of the ultrastrong (USC) and deep strong 
(DSC) coupling regimes of light-matter interaction through the quantum simulation of their dynamics in stan- 
dard circuit QED. The method makes use of a two-tone driving scheme, using state-of-the-art circuit-QED 
technology, and can be easily extended to general cavity-QED setups. We provide examples of USC/DSC 
quantum effects that would be otherwise inaccessible. 
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Introduction. The Jaynes-Cummings model (JCM) fT"! is 
a cornerstone of the field of quantum optics. It describes the 
interaction between a quantized electromagnetic (EM) field 
mode and a two-level atom under two important assumptions. 
First, the interaction is accurately modeled by a dipolar cou- 
pling between the field and the atom. Second, one can apply 
the rotating-wave approximation (RWA) because the coupling 
is small enough when compared to the sum of frequencies of 
the two-level atom and EM field. These restrictions yield a 
solvable model where atom and field exchange one excitation. 

To study experimentally the physics of the JCM, the inter- 
action needs to reach the strong-coupling (SC) regime. This 
can be done by isolating the two-level system from free space 
by means of a cavity with highly reflecting miiTors, making 
the coupling strength much larger than the spontaneous emis- 
sion rate and the cavity decay rate. This type of setup is known 
as cavity quantum electrodynamics (cavity QED) [TTl. Many 
relevant features of the JCM have been successfully tested 
in actual experiments using cavity-QED technology. For in- 
stance, the observation of the vacuum Rabi mode splitting in 
the optical domain with Alkali atoms was reported [4], while 
in the microwave regime, vacuum Rabi oscillations using Ry- 
dberg atoms have also been realized ||3j |5l . 

In 2004, an important step forward was made when 
an artificial two-level atom (or a qubit) was shown to be 
strongly coupled to the EM field inside a superconducting ID 
transmission-line resonator |6|. The newly born circuit QED 
technology was rapidly recognized as a promising architec- 
tural platform from which the study of light-matter interac- 
tion can be extended |7|. Although most circuit QED imple- 
mentations were restricted to the SC regime of the JCM, key 
experiments showing the breakdown of the RWA have been 
recently realized in semiconductor microcavities fE\ and cir- 
cuit QED |9, 10 1 . They have opened up new directions of 
research into the ultrastrong coupling (USC) regime of light- 
matter interaction fTT'T2|, where the RWA can no longer be 
used, leading to novel features such as the creation of pho- 
tons from the quantum vacuum ifTTl . Even though these works 
show that reaching the ultrafast dynamics is feasible, its con- 



trollability becomes very demanding as the light-matter cou- 
pling increases 113]. 

In this work, we introduce the quantum simulation of the 
USC/DSC dynamics in circuit QED with a qubit-cavity sys- 
tem in the SC regime. Our treatment makes use of a novel 
two-tone orthogonal driving to the qubit. We show through 
analytical and numerical calculations that our proposal will 
have access to the regimes of USC (0.1 < g/tu < 1, g/cu be- 
ing the ratio of the coupling over the resonator frequency) and 
DSC lfT4l (g/uj > 1). This will pave the way for the imple- 
mentation of a quantum simulator LISJ for a wide range of 
regimes of light-matter coupling lfT6l in systems where they 
are unattainable. As we will discuss below, this includes the 
simulation of relativistic quantum phenomena, Dicke/spin- 
boson model, Kondo physics, and Jahn-Teller instability 1 17|. 
Although we present our method in the language of circuit 
QED, it can also be realized in microwave cavity QED 1 3 , 1 8 1 . 

The model. The physical system we consider consists of 
a superconducting qubit strongly coupled to a microwave res- 
onator mode. If we work in the qubit degeneracy point, the 
Hamiltonian is given by ||6l 
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where ujq, uj are the qubit and photon frequencies, g stands 
for the coupling strength: Likewise aia^) represent the anni- 
hilation(creation) operators of the photon field mode, whereas 
a, - a+ + a - |e)(.g| + |.9)(e|, o, = |e)(e| - \g){g\, be- 
ing \g) , |e) the ground and excited eigenstates of the qubit. In 
a typical circuit-QED implementation, this Hamiltonian can 
be simplified further by applying the RWA. Precisely 1 19 1, if 
{|a; — ojgl, g} <C w + Wq, then it can be effectively approxi- 
mated as 
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which resembles the celebrated JCM of quantum optics. Per- 
forming the RWA amounts to neglect counter-rotating terms 
a a and cr+a^, leading to a Hamiltonian Eq. ^ where the 
number of excitations is conserved. 
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FIG. 1. Pg{t) obtained by integrating the exact (solid line) Eq. pi and the effective (circles) Hamiltonian Eq. Jm. With the set of parameters 
in the text, fli — 2tv x 700 MHz, we have two cases: (left panel) fl2 = 0; (right panel) Q2 = 2it x 10 MHz. The simulated ratio is 



This Hamiltonian (|2]i will be the starting point of our deriva- 
tions. Consider that the qubit is now orthogonally driven by 
two classical fields. The Hamiltonian of the driven system 
reads 



hio„ 



H = — -(T^ + hwa* a — hg i^a* a + aa*\ 



-hfli e"^iV 



e-^-^iV^l 



hn2 (e"^''a + e"*"^ V^) , (3) 



where fl j , ujj stand for the amplitude and frequency of the 
j— th driving. Note that the orthogonal drivings couple to the 
qubit in a similar fashion as the resonator field does. To obtain 
(J3]l, we have implicitly assumed that the RWA applies not only 
to the qubit-resonator interaction term, but also to both driv- 
ings. Next, we will write ([3| in the reference frame rotating 
with the frequency of the first driving, oji, that is, 



Ti ^ — h— (Tz + h{uj — uii)a* a — hg (cr'a + cra^) 
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1 (a + a^) - hQ.2 (e^^^^-^i^V + g-'^^^-^i^V^) . (4) 



This allows us to transform the original first driving term into 
a time independent one Hq ^ = —Mli [a + cr^^), leaving the 
excitation number unchanged. We will assume this to be 
the most significant term and treat the others perturbatively 
by going into the interaction picture with respect to "Hq ^ , 
n^lt) = e^'^o'^/''' (n^^ -TioA e-''^o't/''\ The physical 

implementation of this transformation based on a Ramsey-like 
pulse is described later in the text. Using the rotated spin basis 
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By tuning the parameters of the external drivings as uji—lu2 = 
2r2i , we can choose the resonant terms in this time-dependent 
Hamiltonian. Then, in case of having a relatively strong first 
driving, Oi, we can approximate the expression above by a 
time-independent effective Hamiltonian as 



-07 ni- ^ t , ^^2 hg 
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Note the resemblance between the original Hamiltonian (fTl) 
and (J6|. While the value of the coupling g is fixed in (J6]), we 
can still tailor the effective parameters by tuning amplitudes 
and frequencies of the drivings. If we can reach values such 
that VL2 ^ (yj — ijJx) '^ g/'i,, the dynamics of the original sys- 
tem will simulate that of a qubit coupled to the resonator with 
a relative interaction strength beyond the SC regime — ideally 
in the USC/DSC. The coupling regime reproduced through 
the effective Hamiltonian (|6]) can be quantified by the ratio 
.gcff/weff, where g^s = 5/2 and ujcs = cj - cji. 

Numerical analysis. To study the feasibility of our pro- 
posal, we have performed numerical calculations with realis- 
tic parameters for circuit QED |6| and compare the evolution 
described by the exact and effective Hamiltonians. In prin- 
ciple, the choice of i7i is unimportant as long as it is strong 
compared to the other parameters involved in this problem. 
Indeed, the evolutions of the Hamiltonians of Eqs. (J3]l and 
(|6]l become more similar with increasing 57i. Nonetheless, 
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FIG. 2. Pg it) obtained by integrating tiie exact Hamiltonian Eq. dSll (liglit solid line, shaded area), and after applying a Ramsey-like pulse 
(dark solid). Both are compared to the effective (circles) Hamiltonian Eq. ^ evolution. The same parameters of Fig. [T] is used: (left panel) 
^^2 = 0; (right panel) ^2 = 27r x 10 MHz. During the Ramsey-like pulse, qubit energy and driving frequency are detuned by — 27r x 200 MHz. 



experimental restrictions are expected to set the limit of this 
driving strength. Throughout the rest of this work, we assume 
fii ^ 27r X 0.7 GHz, which is a realistic value. 

After the discussion regarding the DSC dynamics lfT4l . an 
interesting experiment would be the following. Assume we 
prepare the original SC undriven system in its ground state, 
i.e. |.g,0), and then at time i = we switch on the exter- 
nal drivings. Now the system evolves according to the unitary 
operator which is computed by integrating the driven Hamil- 
tonian Eq. ([3]). For the sake of simplicity, we now assume that 
the evolution of the state is calculated in the rotating reference 
frame of Hamiltonian T-L^{t). Later we will discuss how this 
step can be implemented. The solid line in Fig. [T] shows the 
evolution of Pg{t) for two different cases. Both show a very 
good agreement when compared to the same probability but 
computed from the effective Hamiltonian (circles) of Eq. (|6| 
derived in the strong driving limit. All these simulations are 
done with values of Wg — 27r x 8.01 GHz, cj ~ 27rx 8.01 GHz, 
g = 27r X 20 MHz, wi = 27r x 8 GHz, ^2 = 27r x 6.6 GHz, 
and r^i = 27r X 0.7GHz. Clearly such an amplitude for the 
first strong driving suffices, even to simulate the dynamics of 
a system reaching gcn/^cH = 1. A feature characteristic of 
the DSC dynamics is the existence of a well defined periodic 
evolution for the probability Pg{t), in the case of a degenerate 
qubit. Starting from 1, Pg{t) decays to 0.5 quite rapidly, to 
have a subsequent revival at a time that is equal to the inverse 
of the effective resonator frequency. This is accompanied by 
the generation of photon number wavepackets that oscillate 
in time lfT4ll . Putting O2 = (Fig.[T] left panel) in our sim- 
ulation we observe that Pg{t) presents nearly perfect revivals 
that take place at Wegtrev = ffeffircv = 27r, which corresponds 
to irev = 0.1 /xsec for the set of parameters used. If the ef- 
fective energy of the simulated qubit ^^2 is not zero (Fig. [T] 
right panel) the evolution becomes nonperiodic, producing a 
distortion of the revival peack, which no longer reaches unity. 

Ramsey-like pulse. The computation of the probability 
Pg{t) for the exact Hamiltonian Eq. ^ in Fig. fllhas been 



done in the rotating reference frame used to derive Hamil- 
tonian 7i^{t). However, without going into this interaction 
picture, the evolution of Pg (t) would show a fast oscillating 
term, depicted by the light solid line shown in Fig. [2] Here, 
we propose the following scheme based on a Ramsey-like 
pulse, in order to get rid of this strong oscillation in an ex- 
periment. Imagine that after letting the system evolve with 
the Hamiltonian of Eq. (|3]l for a time t, we switch off non- 
adiabatically both external drivings and apply a detuning of 
about — 27r X 200 MHz to the frequency of the qubit, ujg, from 
its original value 2tt x 8.01 GHz. Next, we switch on a third 
driving with frequency detuned by — 27r x 200 MHz, from the 
value of cji = 27r x 8 GHz, and with amplitude — Oi (op- 
posite phase relative to the first one). The application of this 
Ramsey-like pulse will take the same time t. The wiggly dark 
solid line of Fig. [2] corresponds to the measurement of Pg{t) 
following this method, which matches the evolution obtained 
from the effective (circles) Hamiltonian Eq. (J6]l. 

Nonclassical states and Dirac equation. Nonclassical 
states of the EM field are a paramount resource for quantum 
information processing. However their generation using all- 
optical devices is challenging due to the lack of strong nonlin- 
earities. Through the strong coupling between a qubit and a 
confined EM field, circuit/cavity-QED technology provides a 
way to overcome these limitations ||20l422l . Here we show 
that our setup can be used to generate highly nonclassical 
states of the field mode, without requiring the projective mea- 
surement of the qubit and/or the ability to control accurately 
the qubit-field interaction strength. 

The nonclassicality of a bosonic field can be signaled by the 
Wigner quasi-probability distribution function (WF), defined 
as 



W{a) = lTr(D\a)pfD{a){-lf^), 



(7) 



Pf being the field density matrix and D{a) — exp(Q;a^ — a*a) 
the coherent displacement operator with amplitude a. To 
show the ability of our setup to produce nonclassical states 



of the EM field, we plot in Fig. |3] the WF of the field for 
the same set of parameters studied before, being the evolu- 
tion time set at gcst = tt. The top-left panel on Fig. [3] depicts 
the WF of the state generated when 5I2 = and the qubit is 
measured in its ground state. The result is a well known co- 
herent Schrodinger-cat state with time-dependent amplitude 



(.goff/wcff)(i 



,-iWofft 



1), which is nonclassical. However 



when the qubit is not measured, the state of the field after 
tracing out over the qubit will have as WF the plot of top-right 
panel in Fig. [3] where a classical mixture of coherent states 
with opposite phases is left and any quantum feature is erased. 
When ^2 7^ 0, a surprising property shows up after the qubit 
is traced out: the nonclassicality of the field is not completely 
lost (center-right panel on Fig. |3]l. This effect is more evi- 
dent when Wcff = (bottom-right panel on Fig. [3]). The setup 
might be enhanced by taking advantage of the cavity output 
and novel measurement techniques 1231 . to produce nonclas- 
sical propagating microwaves and lasing in circuit QED ll24l . 
Thecaseof LJcff = Oisphysically most relevant, as it relates 
to the quantum simulation of the 1h-1 Dirac equation |25|. If 
one repeats the derivation from Eq. ([3]l but assumes that both 
external driving come with an additional phase (j) — 7r/2, i.e. 
—hrij (e*("j*+'^V + H.c), and uses the rotated qubit basis 
|±0) = {\g) ± e-"f' |e)) /v^, then HeS becomes 
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for u! — cji, ay = i{a — a^), p = —i{a — a^)/y/2. In the 
appropriate representation, the Schrodinger equation of our 
system is now formally equivalent to that of the 1h-1 Dirac 
equation of a relativistic spin- 1/2 particle. 



in—- — 
dt 



cpay 



• ip, 



(9) 



where p represents the momentum of a Dirac particle of mass 
m and c the speed of light. This shows the ability of our pro- 
posal to access a wide range of physical models. Through 
the engineered Hamiltonian, a qubit at rest coupled to the EM 
field would behave as a massive spin- 1/2 particle moving near 
the speed of light. To observe peculiar effects, such as Zit- 
terbewegung, in our setup one must pay attention at the map- 
ping of the bosonic degree of freedom. While in the original 
Dirac equation (J9]) the operator p corresponds to the momen- 
tum of the Dirac particle, in our Hamiltonian Eq. (jHll this role 
is played by the quadrature p = —i{a — a^)/\/2 of the EM 
field mode in the resonator. Hence, in the simulated Dirac 
equation the dynamics of Zitterbewegung will show up in the 
expectation value of the field quadrature x = {a + a^)/\/2, 
which has the same time evolution as the expectation value of 
the position operator x of a Dirac particle. 

Discussion and conclusions. Although we have disre- 
garded the possible coupling between the orthogonal driving 
and the resonator field, it is easy to show that the effect of 
such a spurious coupling could be sorted out by adding a driv- 
ing to the cavity (acting as a counter coherent displacement). 




FIG. 3. Wigner function W{x,y) of field state after interaction 
time gcff^ = 1", calculated ab initio with parameters of Fig.fT] (Top- 
left panel) 0,2 = 0, Ucs ~ gaH, qubit postselected in ground state. 
(Top-right panel) ^2 ~ 0, lo^h — Qcb, qubit traced out. (Center- 
left panel) JI2 = Ztt x 10 MHz, Uca = ffcff, qubit postselected in 
ground state. (Center-right panel) 0,2 = 2tt x 10 MHz, ojcft = 
Qefc, qubit traced out. (Bottom-left panel) 0,2 = ^tt x 10 MHz, 
i^eB = 0, qubit postselected in ground state. (Bottom-right panel) 
r22 = 27r X 10 MHz, ui^ft = 0, qubit traced out. 



Another source of error stems from qubit dephasing and re- 
laxation rates, as well as the cavity decay rate. Nonetheless, 
interaction times considered in all numerical simulations are 
well below standard decoherence times. To avoid excitation 
of higher levels in the qubit, it is possible to design flux qubits 
where the splitting to the third level is above 30 GHz ll9l [T0l . 
Our method can be extended to the case of a slightly anhar- 
monic qubit via Gaussian shaped DRAG pulses |26| . 



The proposed quantum simulation of a broad range of 
regimes of light-matter coupling may become a building block 
to simulate physics inaccessible in standard quantum optics. 
This includes the Dicke model for multiple qubits, the spin- 
boson model when the resonator is replaced by an open trans- 
mission line, Jahn-Teller instability for several discrete modes 
in the resonator, and relativistic quantum physics such as scat- 
tering of relativistic particles ll27l . Efficient computations of 
these problems may be beyond current numerical capabilities. 
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